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ABSTRACT 

We consider the thermal creation, fusion, evaporation and destruction 
of non-topological solitons (NTS) after a phase transition in the early 
universe. By defining and following NTS statistical equilibrium and de- 
partures from it, we show that depending on particle physics parameters 
one of three possible scenarios occurs. If reaction rates are high enough, 
a period of equilibrium occurs and relic abundances are determined by 
the “freeze-out” temperature. We show that equilibrium first drives most 
NTS’s into their constituents (free <j> particles) and then causes rapid fu- 
sion into large NTS’s. If freeze-out occurs during the first phase, the NTS’s 
are almost entirely destroyed, while if it occurs during the second phase, 
solitosynthesis occurs and NTS’s may be cosmically relevant. For slow 
reaction rates the NTS’s are “born frozen out” and have the abundance 
determined by the phase transition. We develop analytic approximations 
for determining the abundances and test them by numerically integrating 
a reaction network in an expanding universe. Unfortunately, for most of 
the parameter space considered, solito-destruction/evaporation occurs. 
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I. INTRODUCTION 


Non-topological solitons (NTS’s) are classically stable field configurations which 
have been studied by many groups in recent years. 1 They exist in a variety of field 
theories, but most simply one considers a complex scalar field 0 carrying a conserved 
charge coupled to a neutral scalar in such a way to allow regions of false vacuum 
where the 0 is massless and true vacuum where 0 has a finite mass. NTS solutions 
occur when some number Q of the 0’s are trapped in a region of false vacuum and 
are unable to escape because their energy is lower than Qrn the rest energy of Q 
free 0 particles. 

If they exist, non-topological solitons are interesting objects; they could have 
masses ranging from below a proton mass to above a galactic mass, with properties 
quite different from ordinary matter. One is lead to ask whether there is any mech- 
anism for actually forming this kind of coherent state. One possibility, suggested by 
Frieman, Gelmini, Gleiser and Kolb 2 (FGGK), is that during a phase transition in the 
early Universe when regions of false and true vacuum co-exist, a certain number of 0’s 
could be trapped in the false vacuum regions and as these regions shrank they could 
become NTS’s, perhaps surviving until today. FGGK estimated the relic abundance 
of NTS’s, and under several assumptions found that Unts ~ 1 was possible. (Ants 
is the ratio of NTS density to the critical density.) However, FGGK did not consider 
the actual fate of NTS’s after the phase transition. Several possibilities exist; they 
could disassociate or evaporate into free 0 particles, they could absorb free 0’s, fuse 
and become larger, or they could be created thermally by the fusion of free 0’s. The 
object of this paper is to study these mechanisms and the resulting relic abundance 
of NTS’s. 

In many respects the problem of thermal creation and destruction of NTS’s in 
the early Universe is similar to big-bang nucleosynthesis. There, light elements such 
as helium, deuterium, and lithium are synthesised out of protons and neutrons at a 
temperature of around 1 MeV. Borrowing ideas from big-bang nucleosynthesis we will 
take two general and complementary approaches, 1) solving a network of reactions 
involving the annihilation and fusion of a system of free 0’s and NTS’s in an expanding 
Universe, and 2) analytically understanding the results of the network integration by 
defining and following NTS statistical equilibrium and departures from it. 
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By NTS statistical equilibrium (NTSSE) we refer to a state where all reactions 
involving creation and destruction of NTS’s are proceeding faster than the expansion 
rate of the Universe, and the number densities of all species are determined by their 
binding energies and the entropy of the Universe. If NTSSE ever exists, knowledge 
of the abundance of NTS’s created during a phase transition is lost and therefore 
irrelevant. At high enough temperatures and densities we expect NTSSE to obtain, 
but as the temperature drops, the number densities eventually fall so low that reaction 
rates for processes which establish NTSSE become less than the expansion rate and 
the relative abundances of NTS’s “freeze out.” We denote the temperature at which 
this happens as Tp, the freeze-out temperature. Of course, NTS’s can only exist 
after the phase transition finishes, so another important temperature is the Ginzburg 
temperature, Tp, after which false vacuum bubbles are unlikely to spontaneously flip 
into the true vacuum state. If Tg > Tp we expect to have a period of statistical 
equilibrium, while if Tp > Tg we expect the relic abundances to be more or less 
determined by the phase transition as discussed by FGGK. Finally, we note that 
for low enough temperature NTSSE drives all free <f> particles into NTS’s. This is in 
marked contrast with FGGK’s phase transition where for a relic density of NTS’s near 
critical density the relic density of free <f>'s was orders of magnitude larger. We solve 
for the temperature, Tp, at which the density of NTS’s begin to dominate the density 
of free <f >' s and show that if Tp > Tp we truly have solitosynthesis, the creation of 
significant numbers of large NTS’s by the fusion of free <j> particles. For the opposite 
case, Tp > Tp, we find that almost all NTS’s disassociate/evaporate and it is unlikely 
that NTS’s survive in significant numbers. 

We will illustrate these various possibilities using the Lagrangian and phase tran- 
sition model of FGGK, and check the simple conclusions described above by running a 
network of reactions. We find good agreement between the network and the analytical 
approximation for most of the range of parameters we consider. 

The plan of this paper is as follows: In Section II we describe the Lagrangian, 
review some necessary aspects of the NTS solution, and review the phase transition 
scenario of FGGK. In Section III we define and derive the formulas for NTS statistical 
equilibrium. We then develop the important features of NTSSE and find an approx- 
imate formula for Tp. In Section IV we list the reactions which go into maintaining 
NTSSE, write down the network equations and solve them numerically for a trun- 


2 


cated system. We also find an approximate solution of the freeze-out temperature 
Tp and compare this with the numerical results. In Section V we present our results, 
describing conditions under which we have solitosynthesis and conditions under which 
we have solito-destruction. Section VI sums up the paper. 


II. REVIEW OF NON-TOPOLOGICAL 
SOLITONS AND SOLITOGENESIS 

We will use throughout the model and conventions of Frieman, Gelmini, Gleiser 
and Kolb >(FGGK) a . More details concerning non-topological soliton solutions in 
general can be found in Ref. 1. The Lagrangian we consider is 

C = |0^| a + - y(o- J - ao'Y ~ h\<f > | a (<r - <r 0 ) a - y<T 0 (<r - <r 0 ) 3 - A, (1) 

where ^ is a complex scalar field and o’ is a real scalar field. The field <f> has a 
conserved Noether charge Q , a mass of zero at the local minimum <t — <r 0 , and a mass 

= A(o‘_ — o-o ) a at the true minimum <r = <r_. 

By introducing a spherically-symmetric trial solution of the equations of motion 
derived from Eq. (1), and minimizing the resultant energy with respect to the NTS 
radius, one can find the mass (energy) of an NTS of charge Q : Mq = 4tt\/ 2Q 3 ^ 4 A 1 / 4 /3, 
where A is adjusted so that the value of the potential is zero at the true minimum. 
We define the binding energy of an NTS as the difference in mass between the NTS 
of charge Q, and Q massive <f> particles: 

Bq — Qtn^ — Mq. ( 2 ) 

We note that B > 0 (implying that the NTS is classically stable) occurs as long as 
Q > Qnun, where Q m ^ = 1231AA -a (o , _ - o- 0 ) -4 . If Q < Qmm the binding energy is 
negative and at zero temperature the NTS presumably flies apart into free <f> particles. 

For simplicity, we will consider only the value Aj = 0.15Ai used in FGGK, in 
which case the model is described by the parameters Ai, Qmin« and a 0 alone. In this 
case we have A = 0.6A lt r 0 4 , h = 4.24(A 1 /<? mia ) 1 / a , Mq = S-IS^A^V/ 4 , m* = 
5.15onAi a ^ 4 / and the NTS radius Rq = 0.8(Q/\i) l / 4 /<r o . 

In considering the development of the phase transition, FGGK define a critical 
temperature, Tc ~ 2<r 0 , after which the Universe divides up into domains of true 
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and false vacuum separated by domain walls. Because the potential energy density 
of the false vacuum is higher, regions of false vacuum soon shrink, and in this scheme 
the regions of false vacuum which contain a net charge greater than <5 min eventually 
form non-topological solitons. However, not all regions of false vacuum remain false 
vacuum. Thermal fluctuations of the <r field can be large below Tc and regions of 
false vacuum can become true vacuum and vice versa. These fluctuations freeze-out 
at the Ginzburg temperature, To, which FGGK estimate as Tg = 1.30-o/Ai 1 / 3 . They 
estimate the relic density of NTS’s as the density of regions (at T — To) which have 
charge greater than <?min- Defining the number density of free 0’s as n*, the number 
density of free <£’s as and the <j> asymmetry, = (n^ — n^)/n 7 , their estimate of 
the NTS relic number density is 

( \ 3/3 

3-77 — ) expt-QmiaA! 3 /^), ( 3 ) 

^1 min / 


where n 7 = 2 £( 3 )T 3 / 7 r 2 , is the photon number density and is related to the entropy 
density. We see from Eq. (3) that NTS’s become exponentially rare as Qmin increases, 
and that in this scheme almost all NTS have charge Qnu n ; those with Q <5 min are 

suppressed, and of course there are none with Q < Qmin- Equation (3) must break 
down as — ► 0 since the charge in a small region due to Poisson fluctuations will 

then be larger than the charge given by the asymmetry, but we will ignore this here. 3 
The corresponding density of free <j >' s is roughly 

Tty — Tty 


1 - 


CXp( T 'min ) 

3/2 


( 4 ) 


where r m ; n is the factor in the exponential in Eq. (3). Since r m i n must be greater 
than one, 2 as long as > 1 the free <f> density dominates the NTS density. There 
is, therefore, a problem with having a closure density of NTS’s unless the free <f> 
particles (which are massive) can annihilate, or are allowed to decay and the <p's 
inside the NTS’s (which are massless) are not. Even in this case, however, the <f>'s 
inside the NTS can “leak out” quantum mechanically and decay. The calculation of 
NTS decay and the subsequent restrictions on solitogenesis are discussed in Ref. (4) 
and will not be considered further here. This concludes our review of solitogenesis, 
but the interested reader is referred to Refs. (1,2) for more details. 
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III. NTS STATISTICAL EQUILIBRIUM (NTSSE) 


If the system of free <p particles, free <j> particles, and NTS’s is in statistical equilib- 
rium, the number density of each species is determined by its binding energy and the 
temperature and entropy of the system. In kinetic equilibrium the number densities 
of NTS’s of charge Q , and of free cp's and <^’s are given by 

n Q = ^TM’exp^/T^Mg/T) 
n <t> = ^TmJexpO^/T^m^/T) 

n* = exp(/i^/T)A'j(m^/T), (5) 

where T is the temperature, the fi's are chemical potentials, and Ki is a modified 
Bessel function of the second kind and second order. We define statistical equilibrium 
to consist of kinetic equilibrium, defined above, plus chemical equilibrium which, if 
present, means reactions proceed fast enough to ensure relations among the chemical 
potentials. For example, in chemical equilibrium, the reaction <f> + <f> «-> <r + <r proceeds 
fast enough to ensure ^ = 0, and chemical equilibrium of reactions such as 

Q<f> -f X, where X is some state with Q = 0 and is an NTS of charge 

Qi implies = ilq. Other reactions imply {lq = —fiQ. These relations allow us 
to specify the system in terms of known binding energies, the temperature, and one 
chemical potential n = fi#. 

In order to find the number densities we must find /*, which may be done using 
charge conservation. The total charge in a comoving volume, J2 3 , will be conserved; 
Qtot = [ n <t> — n<t> + Q(uq — ng)]i2 3 , where Q mmx is the largest charge NTS un- 

der consideration and in principle is infinity. Since the photon density is proportional 
to the entropy and scales as R 3 we can define the total charge asymmetry, 

/ Qnu \ 

V = n? n*-n*+ £ <?(ng-n< 3 ) , (6) 

\ Q=Qmi* ) 

and remark that r\ is constant as long as entropy is conserved. The value of t) is a free 
parameter and can range from zero (equal numbers of ^’s and <£’s) to 10 -9 (roughly 
the asymmetry of the baryons) to around 1/2 (4>'s as numerous as photons and very 
few 4> f s). 
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For most of our work we will use a non-relativistic expansion of K 2 {x) % e~ lc (ir/2x) 1 ^ 2 , 
although this is not necessary. Taking this limit, we find the number fractions 
Yi = ni/rjn 7 in statistical equilibrium to be 

/if \ 3/2 / r \ 0 

*5* - (4) y ^ 

Y? = (ci»7) _1 (Jy) ' ex Pf(^ “ m *)/ T 1 

YP = y; < ‘exp(-2 M /T), (7) 

where Ci = £(3)y^32 Jir « 3.836, and 

Bq = m t Q - Mq as [l - ‘^j T« (8) 

is the binding energy of an NTS of charge Q. The quantity is the same as lg q 
with Y+ replaced by Y^. The charge conservation condition becomes 

Qnu 

y? - fj* + Y. QlY? - ?«■) = i. (9) 

Q mi. 

which, after specifying 77 and T, can be solved for \l and then all number densities 
can be found. 

In Fig. 1 we display the equilibrium number fractions of <^’s <£’s, and NTS’s for a 
system containing NTS’s from Qmi n = 4 to Q mmx = 5 and their anti-NTS’s. We set 
Ai = 1, ri = .01, and the temperature is divided by Tq so that <To scales out. Figure 2a 
shows the total fraction in NTS’s for the same system but several values of rj. Fig. 2b 
shows the effect of varying the size of the system, parameterized by Qmmx, and Fig. 2c 
shows the NTSSE values of ft for the same parameters values as Fig. 2a. These plots 
show similarities with the corresponding plots of light element abundances during 
nucleosynthesis. At a given temperature, typically one isotope (charge) is favored, 
and as the temperature drops, species with different binding energies takeover. Figs. 1 
and 2 can be understood qualitatively from the charge conservation equation and from 
Eq. (7). Equation (9) can be approximated as l£ q [l — e~ 3M ^ T + (rjTY^^^e 3 ^] = 1, 
where the first term is the contribution from Y£ q , the second from V^ q , and the third 
is the order of magnitude contribution from Ynts- 
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Consider a relatively small value of t]. Then at T = To, fi (Fig. 2c) is small, 
~ Y? » 1 (just about cancelling each other) and Ynts is smaller than either. 
As the temperature drops, F£ q a exp[(/i — m$)/T] drops exponentially (// «C m^), 
as does 5^ q w Y^ q . The exp (B/T) term increases, but the combination of the T^ -1 
and (5^ q )^ -1 dominate at first and so Ynts drops also. As the temperature continues 
to drop, fi approaches m^, and Y^ approaches 1, after which time Y£ q levels out 
and ?; q drops exponentially. (At this point Ynts is insignificant, so Y^ eq cannot drop 
below unity and still satisfy the charge equation.) Since 3^ q is no longer dropping, 
the exponential exp(J5/T) eventually overpowers the T®~ 1 factor and Ynts begins to 
increase exponentially, (at this point Y^ q is insignificant) still having little effect until 
it approaches unity. When Ynts gets close to unity, (l^* 1 )^ -1 must drop to balance 
the exp (B/T) factor, and so Y£ q drops away and we are left with Ynts % 1 and 
negligible amounts of anything else. For larger asymmetry the same evolution occurs, 
except that the drop in Ynts is not as deep and the final rise in Ynts occurs sooner. 

The temperature at which Ints starts to dominate (defined as Tp in the introduc- 
tion) is important in understanding the the evolution of the system, and can be fairly 
well approximated analytically. Since the binding energy of the highest charge NTS is 
most important, in a system with highest allowed charge Q m „, the charge conserva- 
tion equation can be approximated as Y£ q — Y^+Y nts = lj where Ints 
D efining Tp as the temperature when Ints = Y* q = 1/2 and noting that Y£ q (and 
therefore Yg q ) is typically very small at Td, we can solve for the temperature 


Tp(Q m tx) BqJTq /jqv 

For the case we considered numerically, Aj = 1, Qmm = 4, and Qmmx — 5, this becomes 
Td/Tq = — .767/[2.08 + 4ln rj + 6ln(T0/Tc)], which is in good agreement with the 
numerical results of Fig. 2. A plot of our approximate values of Td for various values 
of ij is shown in Fig. 3. 

If we are interested in an infinite system, we note that as Qm^x —* oo, 

r P (oo) m . 

To In ' 


For the values of parameters used in Fig. 2 this gives Td/Tg ~ .4; roughly a factor of 
twenty higher than the Q mmx = 5 case and a factor of two higher than the Q max = 40 
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case. In general 7 p(Q + 1) > Tp(Q ) for all Q, implying that the highest charge 
NTS allowed wants to take up all the free <f>' s. This contrasts with nuclear statistical 
equilibrium where, for example, 4 He dominates before ia C. The difference lies in the 
slightly different pattern of binding energies in the nuclear and soliton cases. 

In conclusion, we see that if statistical equilibrium is maintained until 7 p we 
expect all <^>’s to be absorbed into NTS’s, and in particular into the largest NTS 
possible. If NTSSE continued, eventually all (£’s inside the horizon would be contained 
in one large NTS. In reality, the reactions which maintain NTSSE freeze-out, and the 
actual distribution of NTS’s depends upon how many “reaction times” exist between 
Tp and Tp, and whether in fact, Tp > Tp at all. We now turn to consider departure 
from statistical equilibrium and freeze-out. 


IV. NTS NETWORK AND FREEZE-OUT 

The statistical equilibrium described in the last section is maintained at the micro- 
physical level by the annihilation, creation, evaporation, and fusion of NTS’s and <j> 
particles. So long as these reactions proceed much faster than the expansion rate of 
the Universe, equilibrium can be maintained. The reaction rates are determined by 
cross sections and number densities, and in trying to decide whether or not equilibrium 
exists at a given temperature, one must specify these. If we consider only positive 
charge NTS’s, a set of possible reactions is 

tr + c (a) 

(Qn un~l)4> + X (b) 

*«.,.+> + * (c) (12) 

+ X (d) 

where $g nU , etc. indicates a NTS with charge and X stands for anything 

else (e.g., 7, e + e~, qq). This network would actually extend to + 00), but in 

any numerical integration we must truncate it. We find that many essential features 
are present in even this radically truncated system. In addition to the processes of 
Eq. (12), there will be processes such as $Q mil + $Q mil $(2Qmm) + X , which we 

will leave out since they are beyond the maximum charge of our truncated systems. 
Processes of this type would become important if equilibrium existed long after Tp. 


<f> + 4>_ 

*9«i. + ^ 


etc. 
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One might also question whether the reverse of reactions such as Eq. (12b) can take 
place and we will discuss later the effect of leaving out this reaction. 

In order to follow the evolution of NTS and free <j> number densities one translates 
the processes in Eq. (12), into a set of coupled Boltzmann equations. For example, 
the number density of free (f> particles evolves according to 


= -3 Hnt + (— r„ + f a ) + (Qnd, - l)(r 6 - r b ) - (r c - f c ), (13) 


where the —3 Hn+ is the reduction of due to the expansion of the Universe, the 
are the forward reaction rates of Eq. ( 1 2i ) , and the f,- axe the reverse reaction rates. 
A similar equation can be written for n^, nq, n<? mlB +i, etc. and these are displayed 
in Appendix A. 

For a two body process such as Eq. (12a), r a = (cr a |v|)n^n^, where (<r a |v|) is the 
thermally averaged cross section. Likewise, r b = (o‘fc|u|)n0ng ml> , r e = (<r < ,|t?|)no _,. w < t, 
and Td = (<rrf|Tj|)nQ _,. +in,a. For many-body reactions, the rate is also a matrix element 
averaged over phase space and we can still write f b = (d’fc)(n^)^ mi, " 1 njr , where a b is 
no longer, however, a cross section. 

Equations such as (13) can be simplified by noting that in equilibrium (in a 
non-expanding Universe), r*® — f\ q = 0, that is, in equilibrium the forward and 
backward reactions cancel. 5 Using this detailed balance we have relations such as 
(<7c\ v \) n Q mlm n ^ = (^■c|v|)»Q mU +i^x> w here the “eq” superscript indicates the NTSSE 
abundance. This allows us to replace (<r c |v|) in Eq. (13), and pairs such as r c — r c can 
be written 


r e -f e = (o’ e |t>|) 


n Qm i. n * ~ 


n Q,i. n » n Q-i.+i 


n‘ q 

n Qm U+l 


(14) 


where we have assumed that the X remains in equilibrium throughout the time of 
interest. We can simplify Eq. (13) further by changing variables® from to Yi = 
n i/( n -rV) and from time to x = T/Tg . Then using detailed balance on all reactions, 
Eq. (13) becomes 


dY+ 

dx 


CQ\ 


■<<r.M )(r& - r^r;-) 


- iX^Yq - YpyTY^-'HY?)*—' ') 

-<<r.|«|)(r*r, - TO«/i5,)], 


(15) 
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where 


c = 


3\/5C(3)3GTnpj 


— 0. 014677c nipt, 


(16) 


rn.pi is the Planck mass, and we took the number of degrees of freedom, g, ~ 100. 
There are similar equations for Y$, Yq, etc. and these are displayed in Appendix A. 

Next we need to consider the cross sections appearing in the equations. For the 
<t>+4> -* <r + (T process we can use the interaction contained in the Lagrangian, Eq. (1), 


C e ff ~ 1.06(7 O j Ai® 3 + h\4>\ 3 x 7 


( 17 ) 


where x is the <r field shifted to the true minimum of the potential. Using this we 

find 7 at tree level 

0.014A x 1/a 


(<r«|v|) = 


Q mi n 1 / 2 <r 0 J 


f l - 


A x l ' a <?. 


1/2 


12.5 


(18) 


For processes such as <f> + $q ml , — ► $Q mU +i we are at a loss to calculate the cross 
section, so as a rough approximation, we will set (<r|v|) ~ irR^, where Rq is the 
radius of the NTS taking place in the process. Thus we have (<rfc|v|) ~ (<r c |u|) ~ 
2<? m in 1/a /(Ai 1/ Vo 3 ) and ~ + I/(Ai 1/a cr n a ). Note that all the cross 

sections scale like < 7 0 -a . 

The last step before integrating the coupled set of ODE’s is to specify the initial 
conditions. We will start at Tg in NTSSE, although it would probably be more correct 
to start with the distribution of Ynts [Eq* (3)] derived from the phase transition. 
However, the equations are extremely stiff and we either have a period of statistical 
equilibrium after To, in which case Ynts immediately evolves to YJJjg, or we do not 
have a period of equilibrium, in which case the reactions do nothing and we know 
that we are left with the initial Ynts produced in the phase transition. 

Fig. 4 shows an example of an integration of a network consisting of 0’s, 0’s, and 
NTS’s of charge Q m ^ = 4 through Qm»x = 5. Values of A x = 1, rj = 10 -3 and cr 0 = 10 3 
TeV (4a) or <r 0 = 7 x 10 7 TeV (4b), were used. Note in Fig. 4a that all abundances 
trace their NTSSE values (shown in Fig. 1) until around T/Tq ~ 0.1, when Ynts 
(the sum of the charge four and five NTS abundances) “freezes out” and becomes 
constant. Y^ and Y* follow NTSSE for a good while longer in Fig. 4a, but with the 
smaller cross section (oc <r 0 ~ 3 ) of Fig. 4b they too freeze-out and become constant 
by T/Tg ~ .02. This is the generic picture. Abundances follow their NTSSE values 


10 


until freeze-out, which is determined by the cross sections, after which they become 
constant. 

Since we are potentially interested in systems with large Q max , and since we cannot 
numerically integrate such systems 8 , we would like an analytic approximation for T F 
just as we have for Tr>. Toward this end we first note that 4> and <j> typically stay 
in equilibrium longer than NTS’s so that during freeze-out of the NTS’s we can 
approximate Y+ and by their NTSSE values. The key reaction is then Eq. (12b), 
which is our source of thermal NTS’s, and the equivalent of Eq. (15) for Yq can be 
approximated as 

dV 

-jf * - Yq 1 )- (19) 

Defining A = Yq — >Q q , Eq. (19) can be rewritten as dA/dx = — C7j(<rj,|t;)) Y^ eq A — 
dYq^/dx, and as long as the departure from equilibrium is small we can set dA/dx « 0 
and solve for A: 

Defining freeze-out as A = 1.5Y£ q [see Ref. (6)] and using Eq. (7) for Y^ q and Y^ eq , 
we find 

Tf = (jj/Tg) + 4(A i^/Q^) 

Tg In [(.12Ai 1 ^ 8 Qnun 1 ^ 8 *F 2Tn pi/ <r o) (“ \*F ~ Qp/Tg + 4Ai 3/, *Q / Qmin 1 /* — 

(21) 

where xp = Tp/Tq* For Ai = 1 and Q = = 4, this becomes 

( 22 ) 

where fi must be found from the NTSSE calculation. A plot of /i for several 9 values 
of i) is given in Fig. 2c. 

The approximation, Eq. (21), agrees quite well with the results of the numerical 
integration (typically within 10%-20%) for small tj , but becomes worse as tj becomes 
large (tj > 0.1). For example, with large tj, A x = 1 and = 4 the network gives 
Tf/Tq ^ 0.3 almost independent of the value of <r 0 , while the approximation predicts 
smaller values which vary with <r 0 . A plot of our approximate Tp for several values 
of <r 0 and tj is given in Fig. 5. 


Tf ^ ( /*/?g) -f 2.8 

Tg ~ In [(.14*j/ a m jrf /cr 0 )/(-|* / . - 4/t/T c + 11.3) 


( 20 ) 
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V. RESULTS 


We are now in position to try to decide quantitatively the relic abundance of NTS’s 
as a function of our free parameters 77, cr 0 , Ai, and Qmin- We have discovered that 
there are really just three different possibilities, depending upon the relative order of 
three important events: the end of the phase transition (To), the dominance of the 
NTS’s (Td), and NTS freeze-out (Tjr). 

If Tjr > Tq then the NTS’s are born frozen out, that is they never reach NTSSE. 
NTS’s are neither thermally created, nor are they destroyed, and one is left today 
with roughly the same spectra of NTS’s as was created by the phase transition. This 
was the case tacitly considered by FGGK, and so one expects to have IInts fl$, 
and with their assumptions, it is therefore difficult to have a significant density of 
NTS’s extant today. In terms of the free parameters, we find that if <r<> is greater than 
around 10 13 TeV, then Tf > Tq , where <r 0 is roughly the mass of the free <j> and also 
sets the scale for the (f> + NTS cross section. Keep in mind that the Planck scale is 
1.2 x 10 1# TeV. For large asymmetry Tf/Tq is larger, but even the largest asymmetry 
we considered (77 = 0.5) does not allow Tjr > Tq for <tq < 10 ia TeV . 

If Tq > Tf then, at least for a while after the phase transition, we expect to have 
statistical equilibrium. In this case the number of NTS’s created during the phase 
transition is irrelevant, the final density being determined purely by the temperature 
at which NTSSE ends, along with the asymmetry and NTS binding energies. Typ- 
ically the number fraction of NTS’s is lower than at To, drops exponentially for 
a while, levels off at some very low number fraction, and then rises exponentially 
until it reaches unity, after which time almost all free </>’s disappear. The relation 
between the temperature of the final rise (Td) and freeze-out (Tf) determines the 
NTS relic abundance. If the rise happens before freeze-out (Td > Tf) then we have 
solitosynthesis. Large numbers of NTS’s of various sizes will survive until today and 
the number of free <f> particles will be small. If freeze-out happens first (Tf > Td), the 
number fraction of NTS’s will remain at the very low levels of the NTSSE dip, and 
we have soliton destruction/evaporation. In this case the number fraction of NTS’s 
is insignificant and only free <f>’& and $’s remain today. 

The actual values of Tj? and Td depend upon the parameters of the model: cr 0 , 
Ai, Qmin, the asymmetry, 77. Consistency 10 forces Ai to be near 1 and Qmin to be 
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between roughly 10° and 10 4 . As mentioned, we allow <7 0 to vary up to m ^ = 1.2 x 10 16 
TeV, while we consider 77 between 0 and 1/2. 

In general, as increases both Tp and To decrease. But, since To decreases 

faster, we find that as Qmin increases freeze-out occurs while Ints is still low and 
solito-destruction is the result. Therefore to find solitosynthesis we take Qmin as 
small as possible. Decreasing A x also decreases both Tp and To , this time by roughly 
the same amount. Since there is very little range allowed in the choice of A x we mainly 
consider only A x = 1 . The parameter cq has no effect on To/Tq, but as it decreases, 
Tp decreases, so for solitosynthesis a low value of cr 0 is desirable. (Cross sections 
are proportional to o-o -2 , and a larger cross section will keep things in equilibrium 
longer.) Taking a very low value of <To (and corresponding very large cross section) 
does not help much however, because the dependence of Tp on <To is logarithmic. For 
example, a cross section 10 50 times larger than irR 2 would be needed to have To > Tp 
for Q — = 40 and 77 = 10 -5 . 

The value of the asymmetry is important. For values of 77 near the baryon asym- 
metry (77 ~ 10 -9 ) we always have Tp > To and cannot have a substantial relic 
abundance of NTS’s. As 77 increases both Tp and To increase, but Tp increases to a 
smaller extent. Therefore the best hope for solitosynthesis is when 77 is large. 

Finally we must decide what value of <J max to use. Since To(Q + 1) > To{Q), the 
larger the effective size of the system (Q max )> the larger the area of parameter space 
in which solitosynthesis can occur. With the present analysis we have not been able 
definitively to decide the effective size of the system (we cannot numerically integrate 
a very large system), but we will discuss this question in some detail in Appendix 
B. In Fig. 6 we show the regions of parameter space for which each of the three 
scenarios, solitosynthesis, solito-destruction, and “born frozen-out” takes place. We 
choose favorable (for solitosynthesis) parameter values of = 4 and A x = 1 and 
plot the <ro, 77 plane. The dot-dash line shows the boundary between solito-destruction 
and born frozen-out, destruction occuring to the left of the line. The solid lines show 
the boundary between solitosynthesis and solito-destruction for two values of Q mmx 

(Q = 9 and Q = 44 ). Below these lines solito-destruction/evaporation takes 

place. As discussed in Appendix. B, we favor a boundary with a low value of Qm*x 
as it agrees with an estimate of the allowed region of solitosynthesis made using a 
different method. The dashed line shows the division of parameter space found by 
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requiring enough time to generate an NTS of charge Q mmx via one-body reactions. 
This is a necessary but not sufficient condition for solitosynthesis (See Appendix B). 
Taking the Q mtx = 9 line as the boundary, we see that even for the favorable case 
of Qxnin = 4, Ai = 1 and <r 0 = 1 TeV, we find Tjj > 7> only if t\ > 0.1. This is 
very near the degenerate limit and while such a value for rj is not impossible, it is 
hard to imagine it arising in a natural way. Therefore, we conclude that for most all 
values of our parameters we will not have substantial numbers of NTS’s extant today. 
Typically we either are left with the distribution created by the phase transition or we 
destroy even these, although there is a window of parameter space for which NTS’s 
are naturally produced and could contribute, for example, Ants ~ 1* 

Now we briefly consider what happens if processes such as the reverse of Eq. (12b), 
(Q m -a — 1 )(f> —* + <j> are left out of the network. If NTS’s can only be destroyed 

and not thermally created, the “equilibrium” state has Ynts = 9. Starting from an 
NTS distribution after a phase transition, or from actual statistical equilibrium one 
then finds the “evaporation rate” of NTS’s by integrating such a network. The results 
of this integration are shown in Fig. 7. Comparing Fig. 7 to Fig. 4, we see that Ints 
does drop very quickly to zero (and never rises again), while when thermal creation 
processes are allowed, the drop is more controlled and temporary. In this type of 
scenario, the only way to have any relic NTS’s is to have Tp > Tq, that is, to have 
the NTS’s born frozen out. 

Finally, we mention that while the results presented above are in one sense very 
dependent upon the particular NTS model we used, in another sense they are quite 
model independent. The important temperature scale Td depends primarily on the 
NTS binding energy and asymmetry; the freeze-out temperature depends on these as 
well as the cross sections. So while for different NTS models the precise regions of 
parameter space which give rise to solitosynthesis and solito-destruction will differ, 
we still expect the answers to be given by Eqs. (10) and (21), and to be qualitatively 
the same. 


VI. CONCLUSIONS 

In conclusion, we see that for the model of FGGK there are three generic outcomes, 
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depending on the values of the parameters. If Tp > Tq, NTS’s are born frozen out 
and Ints is determined by the phase transition. If Tq > Tp, a period of 

statistical equilibrium occurs, which erases all knowledge of NTS’s formed during 
the phase transition. If Tp > Tf, which occurs only for fairly extreme values of 
the parameters, Ints ~ 1 and solitosynthesis gives rise to large, perhaps cosmically 
significant abundances of NTS’s. If Tp > Tp, then all NTS’s formed during the phase 
transition are destroyed, Ints *C Y#, and the relic abundance of NTS’s is probably 
insignificant. 
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APPENDIX A 

In this appendix we list the complete set of coupled Boltzmann equations used in 
running the solitosynthesis network. The symbols were all defined in Sections III and 
IV. 
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where <T; stands {or (<7;|u|) and Q = Qnu a throughout. 
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APPENDIX B 

The discussion in the main body of the text which decided the boundary between 
solitosynthesis and solito-destruction/evaporation was very simplistic, but we feel 
probably adequate. In this Appendix we discuss the caveats and our reservations in 
more detail. 

In comparing To with Tp the question arises as to what size system (effective Q mmx ) 
is relevant. It was seen that To(Q + 1) > To{Q), for all Q and therefore the larger 
Qmmx the larger the region of parameter space which would allow solitosynthesis. This 
is basically a bottle-neck question, and since it is exponentially sensitive is difficult 
to answer with confidence. Consider a temperature between To(Qmm + 1000) and 
Td{Q min + 1). NTSSE would drive Ints to unity in a system of size + 1000, 
but force Ynts < 1 in a system of size (?nun + 1. Can the system generate the 
large number of NTS’s needed to reach NTSSE when y eq <j mim , y eq g mim+ i, etc. are 
extremely small and since the $Q m „’s are generated by a sequence of reactions such 

“ b + 4>^ *Q m u+l? 

One approximate way to answer this question is to consider the number of reac- 
tions such as $Q mU + <f> — ► $Q mi .+i which could have taken place between Tq and 
To(Qmmx)‘ This number, N , must be greater than or equal to Q mmx for NTSSE to 
obtain. N is overestimated by IV ~ (crv)n e<l (Q m i*, To(Qmin))tDi where to is the time 
which corresponds to To and n eq ((2min) is evaluated in a system of maximum charge 
Qxnin* (This overestimates N since n C Q mlt would be smaller for a system with a larger 
maximum charge.) Setting N = Q mmx and solving for To(Q max ) one finds the effective 
size (maximum possible charge) of the NTSSE system as a function of the Lagrangian 
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parameters. Q m *x = N > Qmin is a necessary (but not sufficient) condition for soli- 
tosynthesis. In Fig. 6, we show the line <? mai = Qmin (for Qn«in = 4) in the T), (Tq 
plane (dashed line). Note that it falls between the Q m *x = 9 and Q m *x = 44 solid 
lines found previously and argues for a small effective Q mmx . 

This might have been anticipated. At temperatures where Ynts ^ 1, high Q 
NTS’s are very much suppressed compared to NTS’s of charge Q m so one has an 
effective maximum charge of (or Qmin + e, where e is small). If the system freezes 
out before To(Q mmx — Qmin + e), then there was never a time when large numbers of 
NTS’s of charge Q m ^ Qmin + 1, etc. existed, and so no way to generate large numbers 
of Qmax charge NTS’s even if NTSSE would have liked it. 

Another over-simplification of the discussion in the main body of the text was the 
use of Tp, the freeze-out temperature of reaction Eq. (12b), (Qmin— 1)^ * — ♦ +$• 

If any of the reactions, $Q mi ,+Q-i 4 • <j> -+ $q mi ,+< 3 , freeze out then the system will 
stop building. Defining a freeze-out temperature Tp (for example, for the reaction 
above) and using the same method as in Sec. IV. we find 

Tl = (Bq-B Q -i)/T g 

To In [(.057 C<? (<? - l) 1 / J (*?) 1 / J m p i/<r 0 ) (-§*? + Q(-p/T 0 + 4A i 3/4 Q^ 4 ) - B Q /T a )\ ’ 

where Xp = Tp/Tc and 



We find that 7$ > I# +1 for aU Q, so the higher charge reactions freeze-out first and 
are therefore the only ones which need be considered. However, in comparing Tp 
with the previously defined Tp, we find Tp slightly larger or equal to Tp for moderate 
values of Q. Since Tp“ u , Tp“ ta+1 , etc. freeze-out later, the use of Tf, as was done 
in the body of the text, seems adequate. Other limitations of our analysis include 
neglect of the NTS surface energy (clearly not well founded with the small Qmin’s 
considered here) and the use of the zero temperature form of the potential. 

Finally, a potentially serious flaw in our analysis is that the reactions of eq. (12) 
are perhaps not the relevant ones. One might expect many-body reactions such as 
Qmin<i> — ► $q mim + X to exist, as well as $q mU + <j> + <f> — * $Q mi ,+a> etc.. If these go 
at appreciable rates large NTS’s could be built much faster and a larger region of 
parameter space might allow solitosynthesis. 
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10. In order for perturbation theory to be valid (whatever that means in a theory like 
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FIGURE CAPTIONS 


1. Abundances as a function of temperature in NTS statistical equilibrium for 
a system consisting of 4 >, <j > , #4, §s> $4> $5 (Qnm = 5). Parameter values 
Qmin = 4, Ai = 1, and 77 = .01 were chosen. Temperatures are divided by 
Tq so the cro dependence scales out. Since charge is conserved the algebraic 
sum of all the Y ' s (which are actually QT^s) is unity. 

2. Abundances of non-topological solitons in statistical equilibrium for a sys- 
tem with Qmin = 4 and Ai = 1. In (a), the total abundance of NTS’s (Tnts) 
is shown for Qmai = 5 and several values of the asymmetry (77 = 0.5, 10 -2 , 
10 -5 , 10 -9 ). In (b), Ints is shown for 77 = 10 -2 and several values of the 
system size (Q ma* = 5, 40, 200). In (c), the value of the chemical potential, 
ft, is shown for Q m*x = 5 and several values of 77. 

3. The “dominance” temperature, Td, after which NTSSE drives Ynts > 

as a function of the effective size of the system Qmax, and for several values 
of the asymmetry (77 = .5, 10 -2 , 10 -5 , 10~ 9 ). 

4. Abundances of NTS’s found by integrating a truncated network of reactions 
in an expanding universe (to be compared with the equilibrium values dis- 
played in Fig. 1). Parameter values Qmin = 4, Qm.-r = 5, A] = 1, and 
77 = 10 -2 were chosen. Note that anti-NTS’s are not included. Two values 
of <tq (10 3 TeV in (a) and 7 x 10 6 7 TeV in (b)) are shown. The larger cross 
section (oc <to~ 2 ) in (a) allows equilibrium to maintained to lower temper- 
ature, while in (b) freeze-out is seen for Y^, Y^, and Tnts (sum of charge 4 
and 5 NTS abundances). 

5. The “freeze-out” temperature, Tp, defined in the text, as a function of ctq 
for several values of asymmetry (77 = .5, 10 -2 , 10 -5 , 10 -9 ). Parameter 
values Qjoin = 4 and Ai = 1 were chosen and Tp is nearly independent of 
the size of the system, Qm&x* 

6. Parameter space for solitosynthesis, the <to, 77 plane for Qmin = 4 and 

Ai = 1. To the right of the dot-dashed line NTS’s are “born frozen-out” 



and equilibrium never obtains. Below the solid (and dashed) lines solito- 
destruction/evaporation occurs, that is, NTSSE erases knowledge of the 
phase transition, but freeze-out occurs when the abundance of NTS’s is 
approximately zero. Above and to the left of the solid (and dashed) lines 
solitosynthesis occurs. Here, large numbers of NTS’s are synthesized possi- 
bly leading to cosmically relevant abundances. Several methods of deciding 
the boundary between solitosynthesis and solito-destruction are displayed 
(see text and appendix B) but the solid line labelled Qm&x = 9 is probably 
the most relevant. 

7. Pure evaporation, the abundance of NTS’s found by integrating the network 
when “creation” of NTS’s is disallowed. The total abundance of NTS’s 
drops quickly to zero. To be compared with the case when all reactions are 
allowed (dashed line and Fig. 4). 
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